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Abstract 

We develop a method to obtain the large iV renormalization group flows for matrix 
models of 2 dimensional gravity plus branched polymers. This method gives precise 
results for the critical points and exponents for one matrix models. We show that it 
can be generalized to two matrices models and we recover the Ising critical points. 
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1 Introduction: 



Random matrices are useful for a wide range of physical problems. In particular, by the 
means of Feynman rules, random matrices can be interpreted in term of two-dimensional 
surfaces, which themselves are related to two-dimensional quantum gravity 0. 2D quan- 
tum gravity models can be coupled to matter fields, with a non-zero central charge C. 
In term of matrices, this leads to consider multi-matrix models, which, studied near their 
critical points, in the scaling limit, allow to recover a continuous theory. 

Some of the random matrices models have been solved exactly [[J and the continuous 
limit one recovers in the neighborhood of their critical points has been related to C < 1 
quantum gravity models. Although the behaviour of C < 1 models is well understood, 
exact methods have not been able to solve C > 1 models yet. There is a "C = 1 barrier", 
which prevents us from using C > 1 matrix models as a tool to understand C > 1 two- 
dimensional quantum gravity. 

Whereas exact technics are unable to deal with C > 1 models, an approximate method 
is likely to succeed better. In 1992, Brezin and Zinn- Justin [Q] introduced a new method 
to solve random matrices models: the large N renormalization group (RG), where the 
rescaling parameter N is the size of the matrices. Integrating over part of the matrices to 
reduce N one obtains a RG flow in the space of actions. Fixed points should correspond to 
critical points of matrix models in the large N limit and the scaling dimensions of operators 
give the corresponding critical exponents. For instance, the scaling dimension Ao of the 
most relevant operator is related to the "string exponent" 7 S by the relation 7 S = 2 — 

While the large-iV renormalization group method was introduced to study directly 
C > 1 models, its application to already solved (by exact methods) C < 1 models is 
also useful. Indeed, in [|J it was argued by one of the authors that the understanding 
of the behaviour of the flows for C < 1 models would throw light on what happens at 
C > 1: taking into account "branching interactions" in matrix models, it is expected that 
in addition to the gravity fixed point there is another fixed point (corresponding to the 
"branching" transition between 2D gravity and branched polymer behavior). The value 
of the critical exponents when C — > 1 (known from the exact solutions and KPZ scaling) 
support the conjecture of M that at C — 1 these two fixed points would merge, and that 
there would be generically no gravity fixed point for C > 1 models. 

In this paper we develop the large N RG method in order to study precisely ma- 
trix models for 2D gravity + branched polymers. We shall describe first in Sec. ^] the 
renormalization group method as it was introduced by Brezin and Zinn- Justin, and the 
improvements made by Higuchi et al. [H , using the equation of motions. Then we explain 
why these methods are not sufficient to study our class of models. We show in Sec. ^| that 
the linear approximation of [H is not sufficient and we propose a method to go further, 
and we apply it to the one matrix model. The method still requires numerical analysis, 
which is performed in Sec. |], where the resulting RG flow equations are analysed and the 
results of the method are compared to the previous results of || and to the exact results. 
Finally in Sec. |5| we shall consider the generalizations of our method, in particular to the 
Ising plus branched polymer two-matrix model. 



1 



2 The Renormalization Group for Matrix Model 



2.1 General Idea 

We shall consider a matrix model for gravity plus branched polymers. The partition 
function is 

Z = J d<$> e ~ N2v ^ $ hermitean matrix N x N (1) 



where 



V{$) = T 2 + gT A + ^T 2 2 and T n = 1 tr 



(2) 



n 

In Eq. (Q) (7 is the gravity coupling constant and h is the branched polymer coupling 
constant. This model can be solved exactly, thus it should be a good test of renormalization 
group methods. When h is equal to zero, it is the pure gravity model that was considered 
by Brezin and Zinn- Justin in [|J when they first tested their renormalization group (RG) 
method. The general idea of RG is to start from the potential V(<&') for a N' x N' matrix 
$' with N' = N + 1, and to integrate over the last row and the last column of the matrix, 
leading to an effective action V for the N x N remaining matrix $. Performing a linear 
rescaling of $ (wave-function renormalization) to keep the tr($ 2 ) term unchanged one 
obtains the renormalized action V($) + N~ 1 AV($>) and the RG flow equation in the space 
of action n-M^V = AV. 

In it was shown that at order g 3 the action (0) (with h — 0) stayed closed (up to 
a additive shift coresponding to a tr(l) operator), i.e. no new operators were generated 
by the RG. The corresponding beta-function for g was found to have a zero for g < 0, 
corresponding to the critical coupling where pure gravity scaling is recovered. However at 
that order the first numerical result of were only very qualitative, since the error on 
the critical coupling itself is of order 100%. This method was then developped further by 
several authors 0. 

Shortly after, Higuchi, Itoi, Nishigaki and Sakai improved the RG method by using 
the equations of motion of the matrix model (the so-called loop equations) to eliminate 
some of the new operators which appear in the RG transformation at higher orders in g. 
They first considered the pure gravity model with cubic action 

V(1>)=T 2 + gT 3 (3) 

and showed that in the planar limit, the RG equation for the free energy F(g) can be 
written as 

OF OF 

where G is a non-linear function of g, the gravity coupling constant, and of ^ = (T3). 
This is a non linear differential equation, at variance with the standard RG equations which 
should be linear and of the generic form given in [[| 

dF dF 

N M -f3(9)^ + l(g)F = r(g) (5) 
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(where 7(g) and r(g) are regular functions in g). In order to obtain a standard renor- 
malization group equation the authors of || truncated G(g, to the first order in 
For this problem this approximation is in fact quite good. The results they obtained for 
the cubic action @ are : a 0.2% error for the value of the critical coupling and a 2.5% 
error for the eigenvalue corresponding to the string critical exponent 7. Moreover, they 
have generalized this method to a one-matrix model with two couplings: g^T 3 + g^T^ + T 2 
and also to the Ising two-matrix model, with in both cases results quite close to the exact 
results. In particular for the quartic action (g 3 = 0) the relative errors on the coupling and 
the eigenvalue are 1% and 5%. 

Nevertheless, the linear approximation method has only been applied up to now to 
gravity models without branching interactions. The case of gravity plus branched polymers, 
which are the models one has to consider if one wants to verify the scenario of H, would 
be more difficult to treat. One would have to introduce partial derivative of F with 

rjj-, Ull rj J-, 

respect to h, and linearly develop the equations in ^ and ^j. Moreover, the success of 
the linear approximation method of || may be attributed to the fact that the value of the 
critical coupling for pure gravity, g c ~ —0.22, is in fact quite small, since the corrections to 
the linear approximation turn out to be of order g 3 , and are therefore small. For branching 
interactions, we expect that this will not be the case. Already for the pure branched 
polymer model (g = 0, x < 0) the critical point is h c = —0.5, thus we can guess that the 
estimates for the general critical points would be far less precise than those for the pure 
gravity critical point. 

2.2 The Equations of Motion : 

Since in this paper we shall also use the equations of motion to write RG equations for 
matrix models with branching interactions, let us briefly recall the general idea. If one 
starts from a simple action V such as (|2|) (which depends only on the two operators T 2 
and T4) one obtains a variation AV which is a function of all the T 2n for any n. To take 
all these terms into account, one can try to write the action of the renormalization group 
on the most general partition function. This leads to expressions with an infinite number 
of terms which must be truncated in some way. However we know no general and natural 
truncation scheme, and the simplest truncations lead to a very slow convergence for the 
value of critical points, while the estimates for the critical exponents do not converge at all! 
This can be shown explicitely on the simple example of the pure branched polymer model 
with action V($) = T 2 + f T 2 2 (which is in fact a vector model) as discussed in Appendix 
A. 

Nevertheless one must realize that all the operators are not independent. One is free (in 
perturbation theory) to include in the RG transformation a non-linear reparametrization 
of the field variable (i.e. the matrix $) in the partition function, of the form 

$ _> $ + e ($) , e ($) = Ar-i J2M{Tj}]^ k (6) 

k 
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This induces an additional variation in the effective action of the form 

dV de($) 

This variation does not change the physical content of the RG equations, since (SV) = 
where (■ • •) denotes the "vacuum expectation value" (v.e.v. ) Z- 1 J (!$(■■ .) e - N2v W; these 
are the quantum equations of motion of the model, which express the relations between 
the v.e.v. of the traces. This (unphysical) arbitrariness in the RG equations is known as 
the problem of the "redundant operators" . The idea is that one should use the equation 
of motion to simplify the expression of the RG equations and to reduce the number of 
operators and of coupling constants involved in the RG flows. This is what has been done 
in || to reduce the RG equation to the form of Eq. (Q). In Appendix |A] we show how the 
method applies to the simple branched polymer model. 

We are now going to introduce our method, which aims at calculating the renormaliza- 
tion flows of gravity plus branched polymer models. 



3 Our Method : 



3.1 The Action of the Renormalization Group : 

We are now going to explain our method on the example of the one matrix model : V(3>) = 
#T 4 + T 2 + §T| (More general cases as the Ising plus branched polymer model will be 
considered in section |5| and we shall discuss then the generalization of our method). As the 
use of the equations of motion does not allow us to put the renormalized action exactly 
in the same form as the action we would like to study, we are going to consider a slightly 
more general model : 

V($) = if)(T 2 ) + gT 4 (8) 

where if) is of the form : 

oo 

i> = j2 2 ' u ~% T i ( 9 ) 

with hi = 1 and h 2 = h. We are going to show that when the renormalization group acts 
on this model, the renormalized action can be put in the same form, with a renormalized 
if) and a renormalized g, and h\ always equal to one. 
In the following we denote U the derivative of if> 

dip 



U(T 2 



(10) 



We start from $' 



$ v \ 

„ I a (N + 1) x (N + 1) hermitian matrix, $ is a A^ x N 

hermitian matrix, v is a vertical vector, and a a real number. 

v4 / tr$ /2 



Zn+i = J d$' exp 



-(N + i; 



g 



tr$' 
4(A^ + 1) 



+ ip 



2(N + 1) 



fir 
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can be rewritten, to the first order in — , separating the variables $, v, v* and a, as : 



J <i$ dv dv* da exp 



-.V- (V($) + ^AV($) 



with 



Ay = <?T 4 + 2^ - T 2 C/ + <?— + U— + v*(U + ga 2 + ga<$> + + g^- (12) 

Let us introduce first the auxiliary field a and rewrite the part of the integral containing 
the variables v and v* as : 



dvdv* exp 



-N(v*(U + ga 2 + ga$ + g§ 2 + ga)v - g 



a 



(13) 



By integrating over v and v*, we obtain, 

/ = exp[-iV f-(7y + itr [ln(C/ + £« 2 + + £$ 2 + ^) 



(14) 



We then can use a saddle point method and minimize this expression with respect to a . 
We find an implicit equation for the value a s of a : 



— I tr 

N 
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U 



ga s + g& 2 + ga$ 



(15) 



We now have to integrate over a. This can also be done by using the saddle point method. 
The ^ terms disappear thanks to Eq. (|I5]) , and the saddle point a s verifies : 



o 1 

+ C/a fl + — tr 



2ga s + g® 



U + ga 2 s + 0a s $ 







(16) 



Then we can rewrite, where (■ ■ •) means that the average is done over $ 
Zn+i 



(exp 



N (gT i + 2*l,-TzU + g£+g£ + A-g£ 



(17) 



with 



A= 1-1 tr 



ln(*7 + #a 2 + + #$ 2 + #(x s 



and use the factorization property which is valid in the large N limit (AB) = (A){B) + 
0(N~ 2 ) to express the variation of V: 



.ay, 



{N dN } = {9Ta + ^ ~ T * U + 9 , 



+ 9 



+ iv tr 



ln(£/ + #(a s ) 2 + + g<$> 2 + #((x s )) 



(19) 



From the saddle point equation on a s , and the parity of the action which leads to (tr$) = 
we immediately see that (a s ) = is solution of the averaged saddle point equation. 
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Thus : 



at 



with 



(a s > = (-tr(C/ + <7a s + g&y 1 ) 



(21) 



3.2 The Equations of Motion : 



As discussed in Sec. 2.2, one should use the equation of motion and the freedom to add 
redundant operators to the RG transformation to simplify the flows. We now express the 
equations of motion. The change in variables :$—>$ + e$ 2n+1 , n > leads to the 
equation of motion : 

n 

(((ra + 2)gT 2(n+2) + (n + l)UT 2{n+1) - 2^>(n - i) T 2l T 2{n _ i} ) = (22) 



8=0 



Where the T n are defined as in Eq. (0) 



nT n = -tr ($") 



(23) 



We introduce the function f(z) which is the generating function of the (T 2n ), f(z) = 
-^(tr 1 _^ 2 ) (this function / is almost the resolvent of the model). Using Eq. ( p2|) f(z) can 
be rewritten : 



(g + Uz) + J(g + Uzf - Az\g + Uz + 2(yT 2 *) 
/(*) = ( ^ ) 



2z 2 



Noting that 



-./' 



-9 



we immediately find that the solution of the above equation is 

= 2(T a > 

Integrating /, we also have 



(-tr [\n(U + ga s + g® 2 )]} = (\n(U + ga s 



' U+g<r a f(z) — 1 



Finally, denoting 



_9_ 

U 2 



x = 2T 2 U 



(24) 

(25) 

(26) 

(27) 
(28) 
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we obtain 

dV x x 2 
{ ~ N dN = 9Ta + 2 ^ + lnf/ - 2 + ln(1 _ px) + P Y + 
[t5s , z - p - 2z 2 + J(z - p) 2 - Az 2 (z - p - pzx) 
+ l dZ =2^ > (29) 

This expression, if expanded in powers of T 2 , contains a linear term in T 2 : (we recall 

that ip(T 2 ) = Y^jLi^ l ~^hjT 2 3 ). We would like to keep h\ equal to one. This is possible if 
one substracts the coefficient of T 2 , multiplicated by the expression appearing in the first 
equation of motion : 4gT 4 + 2T 2 U — 1. Then we finally obtain a renormalized action exactly 
in the same form as the original action. 

Beyond this point, we shall suppress, for simplicity reasons, the (• • •) in all our expres- 
sions. Indeed, Eq. means that, if we replace in Z N+ i by the right-hand part of 
our equation, the result will be the same. 



4 Numerical analysis : 
4.1 Expansion of the integral : 

We now have obtained a renormalization group equation containing only linear terms in T 4 
and powers of T 2 . So, starting from the action gT 4 + i/j(T 2 ) (with ip{T 2 ) = T 2 + |T 2 2 + ...), 
we can write (at least in principle) the equations for at|j| and : we just have to expand 
the above expression in powers of T 2 . This method, however, leads to very complicated 
expressions with many non-trivial integrals depending on the parameter g. To be able to 
treat our expression easily, it is better to expand it first in powers of p (p is of order g, as 
U — 1 + hT 2 + ... is of order one). The last integral in Eq. ( ^9|) is then also expanded in 
powers of x, and it remains to expand p = — in powers of T 2 , which is easy. 

One should notice that the expansion in p is not trivial : one cannot simply expand 
the integrand and the bounds of integration, as this would lead to divergent expressions. 
In fact, the integral can be expressed as an expansion in p n and p n In p. What one has to 
do to expand the last integral properly is to treat separately the integral between and 
xpa(l — px) -1 , and between xp^(l — px) -1 and p(l — px) -1 . Let us briefly describe these 
operations : the last integral can be expressed as a sum of two integrals : 



1^77 z — p — 2z 2 + J (z — p) 2 — Az 2 (z — p — pzx) 

V — 7 dz = h + I 2 (30) 



/o -2z 3 
The first one given by 



' ~^7p p 2 x 

h = I P x Bi du with z = u and (31) 

' 1 — px 



-2(1 - pxu) 2 + x(l - px) ((1 -u) + ^ (I - u) 2 + Au(l - px) 2 / x) 
Bl = PX -2(l-px«)3 (32) 
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and the second one given by 



f 1 P 
I2 — I B2 da with z = a (33) 

Jx^fp 1 — px 

-2p(l - af + (1 - px) Upx - a) + yj(px - a) 2 + 4pa(l - a) 2 ) 

^ 2 = ^71 A3 ^ 

—2(1 — a) 6 p 

Before expanding the integrand in powers of p, up to the order p 2n , we have to notice 
that, when p tends to zero, the bounds of integration in Ii and I2 tend respectively to 

00 and 0. Let us take I2 as an example. For u ~ large, the terms in the expansion 

of the integrand are of the form p 2n+l u 2n . This implies that the integration leads to a 
term of order p n . Thus, the terms in the expression are really of increasing powers of p, 
which justifies our expansion, keeping in mind that expanding the integrand up to order 
p 2n amounts to expand the integral up to order p n . 

Moreover, when expanding the integrand B\, we obtain expressions of the form : 
P(u)((l — u) 2 + 4u/x)~~ where P(u) is a polynomial. The primitives of such terms 
lead to logarithmic expressions, and we finally obtain a p n and p n In p expansion, that is 
to say, in term of the coupling g, we have an expansion in g n and g n ln(g). We would not 
have realized the existence of these logarithmic terms if we had used a finite number of 
equations of motion; they appear here because, by the use of the f(z) function, we have 
used an infinite number of equations of motion, and expanded properly only after. This 
phenomenon should not come as a surprise, it was already observed in ||. 

Once the expansion in g has been done, we expand U in powers of T2 and, finally, we 
obtain a quite simple expression with two orders of expansion : X in g and Y in T2. 

After development in p at order X = 4 in p, for example, the expression of the integral 

1 = I 1 + I 2 is: 

x - ln(p)^) +p 2 ( - 12 + 8a; - 61n(p) +4xln(p)) + 

/£ ,3/_^l + 84a ._ x 2 + -x 3 - 50ln(p) +36xln(p) -3x 2 ln(p)\ 

4 / 3437 2960 „ 2 4 , 1 , t , . . ' (,>)) 

+p 4 — + -— x - 17 Ax 2 + -x 3 + -x 4 - 490 In p 

V 3 3 3 d 

+400xln(p) - 60x 2 ln(p)) + o(p 4 ) 




Inserting the expansion for / in Eq. ( p9|) we obtain the RG flow equations for the 



couplings g and hj (j < Y). This flow equations can easily be integrated numerically. 



4.2 The shape of the RG flows : 

In figure [l] we show the results obtained at orders X = 5 and Y — 6. The axes represented 
correspond to the two couplings that interest us : the gravity coupling g and the branched 
polymer coupling h = h 2 . Of course, there are Y — 2(= 4) other couplings which are not 
represented there. The RG flows represented here correspond to initial condition hj = 0, for 
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-h 



Figure 1: Renormalization flows, where —h is in abcissa and — g in ordinate. Trajectories 
may seem to bifurcate : in fact, these are two distinct trajectories originating from two 
points near each other but separated by the critical line. 

all the j > 2, i.e. we study the evolution of the model with inital action = gT 4 + T 2 + |T| 
as in Eq. (|2|). It is easy to find the critical line in the (g, h) plane: it separates the domain 
where one flows towards the Gaussian fixed point (g — > 0, h — > 0) from the domain where 
g and h diverge. We have chosen here to show only the flows for (g, h) near the critical 
line, on both sides of it. Since under the RG flows the hj's becomes non-zero, and since 
we project the flows on the (g, h) plane the flows may seem to cross, of course this is 
unphysical and a projection of what happens in higher dimension. 

First one recovers the correct phase diagram and critical line (with an average 5% 
relative error at this order). The critical line is separated into two parts by a multicritical 
point C. Starting from the rightmost part of the critical line the flows are driven towards 
a fixed point B' which lies on the g = plane. Starting from the leftmost part of the 
critical line the flows are driven towards a fixed point A' (with g, h, h^, ... ^ 0). We 
have represented on Fig |I]the "pull-back" A of A', i.e. the point on the critical line which 
flows in the fastest way towards A' (this is equivalent to say that the leading subdominant 
corrections to scaling vanish at A). Both fixed points B' and A' have one unstable direction. 
Thus B' corresponds to branched polymers and A' is the fixed point corresponding to 
pure gravity. Starting from the point C one flows towards a fixed point C (also with 
with g, h, hz, ■ ■ ■ ^ 0) with two unstable directions, which should therefore correspond to 
the "branching transition" between pure gravity and branched polymer scaling behaviour. 
These features of the RG flows are in full agreement with the picture which was proposed 
in ||]. This agreement is confirmed by studying more precisely the numerical results. 

4.3 Quantitative results : 

At this order of approximation (X = 5, Y = 5), we obtain for the critical point g c located 
on the axis fe = 0a precision which is smaller than the one obtained by Higuchi et al.. Our 
computations, however, can be performed for several values of X and Y, and the results 
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can be extrapolated to obtain a very good precision. For example, for g c , critical value of 
g for h = 0, we have : 



X 


2 


3 


4 


5 


oo 


9c 


-.095415(5) 


-.090745(5) 


-.088725(5) 


-.087585(5) 


-.083325(5) 



.094 
. 092 



.088 
0.086 
. 084 




Figure 2: Extrapolating the h = critical point — g c as a function of 1/X 



We then extrapolate these results by a polynom of degree three in X 1 . Fig |4.3j shows 



that this extrapolation behaves almost as a straight line : the value of g c converges as 
The extrapolated value for X = oo is —0.083325(5), while the exact value of g c is 
— i ~ —0.083333(5). We have thus obtained a 0.016% relative error, to be compared with 
the 0.9% relative error of |J. 

By the same method, we can obtain the whole critical line in the plane (g,h), which 
was not obtained in ref.||. On the pure branched polymer line (g = 0), we obtain for the 
branched polymer critical point h c ~ —0.50037 instead of the exact value h c = —0.5 i.e. a 
0.07% relative error. For the position of the multicritical point C we obtain respectively 
for h c and g c a 0.02% and 0.5% relative errors. 

Linearizing the RG flow equations in the vicinity of the fixed points, we obtain also good 
results for the critical exponents. For instance, studying the multicritical fixed point C, 
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we obtain for the largest eigenvalue Ao — 1.185 , i.e. a 1.25% relative error when compared 
to the exact result 6/5 = 1.2 ! 

It is not surprising that the critical exponents converge less rapidly than the positions 
of the critical points. Better but more tedious calculations could be done, for we can go in 
principle to any orders X and Y, and the limit X — > oo and Y —>■ oo should give the exact 
result, but this is not our purpose here. 



5 Multimatrix Models 

Our method can be generalized to a two-matrix model which describes the Ising plus 
branched polymer model : 

tr(A 3 + B 3 ) tr(A 2 + B 2 ) h ( ti(A 2 + B 2 )\ 2 ti AB . , 

Vo = 9 3N + 2N + 2{ 2N ) + C °^T (36) 

A and B are two N x N hermitian matrices, g is the gravity coupling constant, h the 
branched polymer coupling constant and c corresponds to the temperature of the Ising 
model. For simplicity reasons, we have not introduced a external field, which would have 
broken the symmetry in A and B. 

As previously, we are led to study a more general model : if we note $=^^ £>^)'^ 1 = 
^ ^ J ^ , and So = Id = ^ J J ^ , the model we are going to study can be expressed as : 

y 3iV * \ 2N 7 2N J K } 

Where ip is, as in the case of one matrix models, a function of the quadratic terms in 
Tf r _ tr(f") _ tr(^"+B") , rp — tr$(ft$fr) _ tr(AB) , a± _ n anH _a±_ _ 

Here c is a series in Tab an d ^2, with constant coefficient Cq. The renormalization group 
equation reads : 

ot r 3 2 

-N— = gT 3 + 2i>- T 2 U - cT AB + <?y + (U + 2c) y + tr ln((C7 + ga)5 + cS 1 + g$) (38) 
Where a — T\. 

To simplify this expression, we have to use the equations of motion. Similarly as what 
was done in 0, we can express, by introducing $ — > $ + e<5 in . . . ^^(zSo + ~8\ + 
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changes in variables, Wq = {jrtr zS + -Si + $ ) as the solution of a quartic polynomial : 



- (* + ~ g )(U - gz)j \-±{U + 2c - gz)W> + (2^ + 
(c + - gz)W + c(U - gz)2T 2 + gcT AB \ - ( ^ - (z([7 - <?z) 



+ — [/ + ^ Wo - c + - (U + 2c - gz c + gz)(U - gz) 

2# g 
(w2 c ^ 

[-Y-V- 9 Z )( Z + -Wo + 2gT 2 + 2(U - gz) 







(39) 



We have then to integrate : 



\Z 

a \ 



Without going through all the details of the procedure, let us note that, to expand W , 
we have to be cautious. Indeed, as for the one-matrix model, the integral will have to be 
cut into two parts, as U — gz can be large [z ~ oo) or small (z ~ ~). This will lead to two 
different changes of variables in the integral and, once more, to logarithmic terms in g and 
in -co. 

Let us also stress that one has to expand the integral both in powers of g and of Co. 
The easiest way to do so is to expand simultaneously / in powers of g and Co, with g 2 ~ c . 
Then, we can compute approximate flows for this model. 

Figure 3 shows approximate RG flows, at x = 0, at order X = 6 in g (i.e. three in Co) 
and Y = 2 in T 2 , where g is in abcissa and Cq in ordinate. 




Figure 3: Approximation of the flows for the Ising model, g is in abcissa, c in ordinate 
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Here, we recover the pure gravity (c = 0, x — 0) critical point at g c ~ 0.262, i.e. with 
19.5% error (the exact critical point is 432 -1 / 4 ~ 0.219). By extrapolating the two first 
results at X = 2 and X = 3, however, we obtain an extrapolated g c ~ 0.222 that is to say 
only 1.2% error! Moreover, we can see on figure 3 that we recover the Ising bicritical point 
C at c ~ —0.14 (the exact value is 0.15), and the shape of the critical line. 

All these results show we can compute good approximations of the flows, not only for 
one-matrix models, but also for multimatrix models. We can theoretically compute that 
way the flows of an open chain of k matrices with nearest neighbour coupling, for k = 1, 
2 (Ising), 3 . . . This series of models is all the more interesting as we know that when k — > oo, 
the central charge of the model C — > 1. These models could thus allow us to verify the 
evolution of the flows with C and their shape when C is equal to one. 

This leads us, however, to a practical problem : though the initial model (the open 
chain) has the same coupling constants for all the matrices of the chain, the action of the 
renormalization group leads to almost as many coupling constants as matrices. This comes 
from the fact the roles of the matrices of the open chain are not symmetric. Thus, we do 
not know if it is practically manageable to study, for example, a k — 10 matrices open 
chain, or if it leads to too long computations. A solution is to study a symmetric problem : 
the closed chain, for example (we do not know its exact solution), or the k- matrices Potts 
model, where all matrices are coupled. The study of the latter model could be indeed very 
interesting : when k > 4, then C > 1 and we would thus enter the C > 1 domain. 

6 Conclusion 

In this paper we have developped the large N renormalisation group method to study ma- 
trix models containing interaction terms corresponding to branched polymer interactions. 
We have shown the analytical basis and the successes of our method. It can deal with 
models containing branched polymers, it gives us the shape of the flows, and also good ap- 
proximations of the position of the critical points and critical exponents of the models. We 
have applied our method to the case of the pure gravity plus branched polymer one-matrix 
model, and to the case of the Ising two-matrix model. Our method is an approximation 
method : the exact expressions must be truncated to a certain order to be numerically 
manageable (the ideal case of the infinite chain being the exact solution). However, the 
extrapolation of the first orders gives good results without taking high computation times. 
But, when studying models with a growing complexity (/c-matrices open chains), we may 
reach for big k the practical limits of the method, and thus it would be a good thing to find 
more technical simplifications. We also plan to study A;-matrices Potts models, which are 
very symmetric models (so technically simpler) but which would allow to cross the C > 1 
barrier for k > 4, and are thus theoretically more complex models. 
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A The Vector Model 



The large N RG method has already been applied to the (very) simple case of the vector 
model in ||. This model corresponds to the limit g = of the matrix model with action 
(§), since then the U(iV) matrix model becomes a 0(iV 2 ) vector model). The authors of || 
found that if one does not use the equations of motion (i.e. if no non-linear reparametriza- 
tion of the fields is performed) the RG flow equation seems to lead to strange results, with 
the apparent existence of a one-parameter family of fixed points. However they showed 
that when using the equations of motion one can reduce the RG flows to much simpler 
flows in a finite dimensional space of couplings. Moreover in that case the flows can be 
easily calculated exactly and it is found that one recovers the exact critical points and 
critical exponents. 

The purpose of this appendix is just to present the results of a simple exercise: how do 
the results of the RG method change if one does not use the equations of motion, or if one 
uses the equations of motion in a approximate way? 

We start from the action (H) with g = for a, N x N Hermitean matrix 

m)=Y,^kTj , T*=±T±{^ (40) 

Since we can rewrite T 2 = |$| 2 where <3> is the iV 2 -dimensional real vector whose components 
are the real and imaginary part of the matrix elements of $, the model reduces to the 
vector model studied in ||. The RG flow equation for the potential ip can be written 
exactly, either by integrating explicitely over some components of $, or by computing 
exactly the effective potential at large N and studying its variation with N. Both methods 
yield, to the first order in 

= 2 ^ ~ 2T ^' + (41) 

This simple RG equation becomes even simpler if we consider instead of ip its derivative 
U = ip' an d if we invert the function U(T 2 ) and consider instead the function X = U^ 1 
defined as T 2 = X{U) (this is at least possible for T 2 small). The RG flow equation becomes 
linear for the function X 

d 1 

N—X = 2X (42) 

dN U K J 

and is trivial to solve. If we start at Nq from the initial potential ipo, i.e. from the initial 

function X with X = Uq~ x = i^'q~ 1 we get at iV = AiVo 

X(U;N) = l -^ + A 2 X {U) (43) 

and inverting again T 2 = X(U) we obtain the renormalized derivative of the potential 
[/(T 2 , N) = -j^ripiTii N). Performing a linear rescaling $ — > Z~ 1 / 2 & to keep the coefficient 
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hi of the T2 term fixed amounts to changing Xq(U) — » ZXq(ZU) in Eq. (|43|), with Z = 
Z(A) a linear rescaling factor fine tuned such that the constraint ip'(0) = hi (i.e. X{h{) = 
0) is kept for N ^ 1. 




RG flows of t7 = V' 

ELS cL function of 
/i. The initial function Uq is the thick 
black line, the large iV limit is the 
thick dashed line. For A > 1 U devel- 
ops a singularity which moves with A 
along the thin black curve. The renor- 
malized U are depicted by the thin 
(color) curves, from the UV with A 
small (blue) to the IR with A large 
(red), for h < h c the singularity goes 
to 00 and U — > 1 goes to the Gaussian 
x fixed point. For h = h c = 1/2 U goes to 
the non-trivial fixed point. For h > h c 
the singularity hits the origin and the 
RG flows diverge in a finite time. 

Figure 4: 

First let us study the exact flows if one starts from the quartic potential ip = T 2 — 1^ 2 2 , 
i.e. from the initial function X = (1 — U)/h. For A > 1 small, the function U remains 
analytic around the origin, but it develops a square root singularity at a finite X = X S (A), 
which starts for A = 1 at the zero X = h~ l of Uq, i.e. at the critical point of ipQ. Of course 
the flow can be studied analytically but they are better depicted graphically (see Fig. |) 

• If < h < h c — 1/2, the singularity X S (N) goes to infinity as — > 00 and the 
function U(X, N) tends towards the constant function U(X) = 1, which corresponds 
to the Gaussian potential ip = T2 (Gaussian fixed point) (Fig. |^.a). 
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• At the critical value h c = 1/2 the singularity X s tends toward a, finite value X* = 1/2, 
so that the function U tends toward a non-analytic fixed point (Fig. f|.b) 

U *( X ) = { undefined if X > X* ^ 

• For h > h c = 1/2, the singularity X S (N) reaches the origin X = in a finite RG 
time A s = 2h/(2h — 1), the potential ip becomes singular at the origin and the RG 
flow thus diverges in a finite time and reaches no fixed point (Fig. |].c). 

This analysis can be extended to general initial potentials. We thus recover (without using 
the equations of motion) a sensible picture of the RG flows: the attraction domain of the 
Gaussian fixed point, (h < h c ) is separated from the domain where RG flows diverge by 
a critical (unstable) manifold where one is driven towards the non-trivial fixed point U* 
(which corresponds to branched polymers). The only subtle point is that the non trivial 
fixed point is non-analytic and cannot be distinguished by a local analysis around the origin 
(T2 = 0) from the Gaussian fixed point. This explains the apparent paradoxes of ||. 

Let us mention that if instead of using the potential U(X) one uses the inverse function 
Y(U) = X(U) - 1/2U, the RG flot is very simple: Y{U) -> A 2 Y(U). The structure of the 
flow is described by the expansion of the function Y(U) around the largest zero U* such 
that Y(U*) = 0, and the (only) relevant scaling field to corresponds to the first derivative 
Y'(U*) of Y (since the critical manifold corresponds to t = 0). This to scales with A as A 2 , 
therefore has scaling dimension 2. However the mapping U(X) — > to is highly non-linear, 
and becomes singular along the critical manifold to — 0. Therefore this does not contradict 
the fact that the real scaling dimension of the relevant operator is 4/3. 



One can try to truncate the RG flow equation Eq. (41) in the most naive way: we keep 
only the couplings hj with dimension j < K in ip, then expand the In [?/>'] in powers of T2 
and truncate this expansion at order K in Ti. We thus obtain for fixed K approximate 
RG flow equations for the K couplings hj, which are of the standard form N-^hj = /3j[h], 
with the P functions polynomials of order K in the h^s. At a given truncation order K 
the explicit form of these functions is not especially illuminating and will not be given 
here. Using computing software the approximate fixed points can be found exactly and 
the structure of the RG flows studied. We find indeed a non-trivial fixed point ip* with 
one unstable direction, which should correspond to the branched polymer fixed point. 

The derivative U* = ip*' is depicted on Fig. | as a function of the order of truncation 
K (U is a polynomial of degree K — 1. One sees that as K increase the approximated 
fixed points converge towards the exact (but singular) fixed point. However a more precise 
analysis (that we do not reproduce here) shows that the convergence is very slow, typically 
as 1/K. The finite K estimates for the critical coupling h c (if all higher order couplings 
hj, j > 2 are set to zero) converge towards the exact value 1/2 at the same rate. 

This very slow convergence is insufficient to obtain good estimates for the critical ex- 
ponents. It turns out that within this approximation scheme, the scaling dimension A^ ppr 
of the scaling fields at the approximate fixed point are independent of the trucation order 
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Figure 5: 

K\ Indeed they are found to be integers (1, —1, —2, —3, . . .). The +1 dimension should 
correspond to the dimension A of the relevant perturbation, which is known to be 4/3. 
Thus, although the estimates for the critical points converge towards the correct result as 
K — > oo, the estimates for the scaling exponents do not! A procedure to accelerate the 
convergence is required. This is precisely what the equation of motions are doing. 
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